Physical Layer Network Coding for Two-Way 
Relaying with QAM and Latin Squares 



Vishnu Namboodiri 

Dept. of ECE, Indian Institute of Science 

Bangalore 560012, India 

Email: vishnukk@ece.iisc.ernet.in 



B. Sundar Raj an 

Dept. of ECE, Indian Institute of Science, 

Bangalore 560012, India 

Email: bsrajan@ece.iisc.ernet.in 



o 

(N 



in 



C/3 



> 

&\ 
^D 
(N 

en 

o 

(N 



X 



Abstract — The design of modulation schemes for the physi- 
cal layer network-coded two way relaying scenario has been 
extensively studied recently with the protocol which employs two 
phases: Multiple access (MA) Phase and Broadcast (BC) Phase. 
It was observed by Koike-Akino et al. that adaptively changing 
the network coding map used at the relay according to the 
channel conditions greatly reduces the impact of multiple access 
interference which occurs at the relay during the MA Phase 
and all these network coding maps should satisfy a requirement 
called the exclusive law. In |10| it is shown that every network 
coding map that satisfies the exclusive law is representable by 
a Latin Square and conversely, and this relationship can be 
used to get the network coding maps satisfying the exclusive 
law. But, only the scenario in which the end nodes use Af-PSK 
signal sets (where M is of the form 2 A , A being any positive 
integer) is extensively studied in |10|. In this paper, we address 
the case in which the end nodes use M-QAM signal sets (where 
M is of the form 2 2A , A being any positive integer). In a fading 
scenario, for certain channel conditions •ye? , termed singular 
fade states, the MA phase performance is greatly reduced. We 
show that the square QAM signal sets give lesser number of 
singular fade states compared to PSK signal sets. Because of 
this, the complexity at the relay is enormously reduced. Moreover, 
lesser number of overhead bits are required in the BC phase. The 
fade state ye j9 = 1 is singular for all constellations of arbitrary 
size including PSK and QAM. For arbitrary PSK constellation 
it is well known that the Latin Square obtained by bit-wise XOR 
mapping removes this singularity. We show that XOR mapping 
fails to remove this singularity for QAM of size more greater than 
4 and show that a doubly block circulant Latin Square removes 
this singularity. Simulation results are presented to show the 
superiority of QAM over PSK. 

I. Preliminaries and Background 

We consider the two-way wireless relaying scenario shown 
in Figfl] where bi-directional data transfer takes place between 
the nodes A and B with the help of the relay R. It is 
assumed that all the three nodes operate in half-duplex mode, 
i.e., they cannot transmit and receive simultaneously in the 
same frequency band. The relaying protocol consists of the 
following two phases: the multiple access (MA) phase, during 
which A and B simultaneously transmit to R using identical 
square M-QAM signal sets and the broadcast (BC) phase 
during which R transmits to A and B using possibly with 
another square M-QAM or constellations of size more than 
M. Network coding is employed at R in such a way that A 
(B) can decode the message of B (A), given that A (B) knows 
its own message. 
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(b) BC Phase 
Fig. 1. The Two Way Relay Channel 



A. Background 

The concept of physical layer network coding has attracted a 
lot of attention in recent times. The idea of physical layer net- 
work coding for the two way relay channel was first introduced 
in HI, where the multiple access interference occurring at the 
relay was exploited so that the communication between the 
end nodes can be done using a two stage protocol. Information 
theoretic studies for the physical layer network coding scenario 
were reported in [2|, [3|. The design principles governing the 
choice of modulation schemes to be used at the nodes for 
uncoded transmission were studied in ||4). An extension for 
the case when the nodes use convolutional codes was done 
in O. A multi-level coding scheme for the two-way relaying 
scenario was proposed in [6|. 

It was observed in flU that for uncoded transmission, the 
network coding map used at the relay needs to be changed 
adaptively according to the channel fade coefficients, in order 
to minimize the impact of the multiple access interference. 
The proposed Latin Square scheme was studied in [10|, ifTTl 
by considering a two way relaying using A/-PSK signal sets 
at the end nodes. In (4) analysis of 16-QAM is done under 
the assumption that precoding is done at the end nodes. We 
address the situation where no precoding assumption is made 
and to the best of our knowledge no work has been reported 



for such a scenario with general M-QAM modulation. 

B. Signal Model 

Multiple Access (MA) Phase: Let S denote the square M- 
QAM constellation used at A and B, where M = 2 2A , A 
being a positive integer. Assume that A (B) wants to transmit 
an 2A-bit binary tuple to B (A). Let fi : S — > F 2 2a denote the 
mapping from complex symbols to bits used at A and B. Let 
/j,(xa) = sa, h(xb) = sb £ S denote the complex symbols 
transmitted by A and B respectively, where sa,sb £ F 2 2a. 
The received signal at R is given by, 

Y R = H A x A + H B x B + Z R , 

where Ha and Hb are the fading coefficients associated with 
the A-R and B-R links respectively. The additive noise Z R 
is assumed to be CJ\f(0,a 2 ), where CJ\f(0,a 2 ) denotes the 
circularly symmetric complex Gaussian random variable with 
variance a 2 . We assume a block fading scenario, with the ratio 
Hb/Ha denoted as z — je^ , where 7 £ K + and — 7r < 9 < 
7r, is referred as the fade state and for simplicity, also denoted 
by ( 7 ,0). 

Let S R ("f,9) denote the effective constellation at the relay 
during the MA Phase, i.e., 

Sfl(7, 0) = {xi + "/e^Xjlxi, Xj £ S} , 

and rf m j n (7e J ) denote the minimum distance between the 
points in S R (j,9), i.e., 



Broadcast (BC) Phase: Depending on the value of 7e J , R 
chooses a map .A/f 7 ' 9 : S 2 — >• S', where S' is the signal set (of 
size between M and M 2 ) used by R during BC phase. The 
elements in S 2 which are mapped on to the same complex 
number in S' by the map M. 1 ^ are said to form a cluster. 
Let {£1, £2, •••) Ci} denote the set of all such clusters. The 
formation of clusters is called clustering, and denoted by C 7e 
to indicate that it is a function of 'ye^ 6 . The received signals 
at A and B during the BC phase are respectively given by, 



Y A = H' A X R + Z A ,Y B = H' b X r + Z B , 



(3) 



where X R — .M 7 ' (a:^, Xb) £ S' is the complex number 
transmitted by R. The fading coefficients corresponding to the 
R-A and R-B links are denoted by H' A and H' B respectively 
and the additive noises Za and Zb are CAT(0, a 2 ). 

In order to ensure that A (B) is able to decode B's (A's) 
message, the clustering C should satisfy the exclusive law |0J, 



i.e., 



Mri' e {x A ,x B )^M"l> e {x' A ,x B ), iorx A ^x' A ,Vx B £ S, 
M"*' e (x A ,x B )^M~t' 6 (x A ,x' B ), iorx B £ x' B ,Vx A G S. 



(4) 



Definition 2: The cluster distance between a pair of clusters 
Li and Cj is the minimum among all the distances calculated 
between the points xa + je^XB^x^ + 'fe^x'g £ Sr(7, 9) 
where (xa,%b) £ £i and (x' A ,x' B ) £ Cj. The minimum 
cluster distance of the clustering C is the minimum among 
all the cluster distances, i.e., 



, (-ye J ) = min | (x A - x' A ) + "/e J (x B 

' )es 2 



(x A ,x B ) ,(x' A ,x' B )€i 
(x A ,x B )^(x' A ,x' B ) 



e 'b)\- (!) d C ■ hei°) 



From ([TJ, it is clear that there exists values of 7e- ? for 
which dminil^ 6 ) = 0. Let H = {^ e £ C|d mm ( 7 ,0) = 0}. 
The elements of H are said to be the singular fade states. 
Singular fade states can also be defined as 

Definition 1: A fade state ^e^ 6 is said to be a singular fade 
state, if the cardinality of the signal set S R 1(7, 9) is less than 
M 2 . 

For example, consider the case when symmetric 4-QAM 
signal set used at the nodes A and B, i.e., S = {{±l±j)/y/2}. 
For jeJ 9 = (1 +i)/2, d min (7e j8 ) = 0, since, 
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Alternatively, when 7e j6 ' = (1 + j)/2, the constellation 
S R (j,9) has only 12 «16) points. Hence je l9 = (l+j)/2 
is a singular fade state for the case when 4-QAM signal set is 
used at A and B. Let (xa,xb) £ S 2 denote the Maximum 
Likelihood (ML) estimate of (xa,xb) at R based on the 
received complex number Y R , i.e., 

(x A , x B ) = arg min \Y R - H A x' A - H B x' B \. (2) 
(x' A ,x' B )es 2 



= min 

(x A ,x B ),(x' A ,x B ) 
65 2 , 
M~<- < *(x A ,x B )j:M~<- e {x' J 



The minimum cluster distance determines the performance 
during the MA phase of relaying. The performance during the 
BC phase is determined by the minimum distance of the signal 
set S' . For values of je^ e in the neighborhood of the singular 
fade states, the value of d m i n (^je^ B ) is greatly reduced, a 
phenomenon referred as distance shortening. To avoid distance 
shortening, for each singular fade state, a clustering needs to be 
chosen such that the minimum cluster distance at the singular 
fade state is non-zero and is also maximized. 

A clustering C is said to remove a singular fade state h £ W, 
if d^„(/i) > 0. For a singular fade state h £ H, let C^ 
denote a clustering which removes the singular fade state 
h (if there are multiple clusterings which remove the same 
singular fade state h, consider a clustering which maximizes 
the minimum cluster distance). Let C-u — {C^h} '■ h £ T~L] 
denote the set of all such clusterings. Let d m i n (C^ h \y : 9') 
be defined as, 



m(C«, 7 ',«')= min I {x A - x' A ) + 7V' {x B - x' B ) \. 

(x A ,x B ),(x' A ,x' B )es , 

M^ h Ux A ,x B )^Mi h 'i(x l A ,x' B ) 



The quantity d m i n (C^ ,7,' 9') is referred to as the mini- 
mum cluster distance of the clustering C^ evaluated at j'e^ 8 . 

In practice, the channel fade state need not be a singular 
fade state. In such a scenario, among all the clusterings which 
remove the singular fade states, the one which maximizes 
the minimum cluster distance is chosen. In other words, for 
"f'e^ 8 $■ H, the clustering C 7 ' 6 is chosen to be C^ h \ which 
satisfies d mm {C {h} ,i ,0') > d min {C^ h '\i ,8'), Vh £ h' € 
W. Since the clusterings which remove the singular fade states 
are known to all the three nodes and are finite in number, the 
clustering used for a particular realization of the fade state can 
be indicated by R to A and B using overhead bits. 

Example 1: In the case of BPSK, if channel condition is 
7 = 1 and 9 = the distance between the pairs (0, 1)(1,0) 
is zero as in Fig|2a).The following clustering remove this 
singular fade state. 

{{(0,1)(1,0)},{(1,1)(0,0)}} 

The minimum cluster distance is non zero in this clustering. 
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(b) Effective Constellation at the relay for -y = 1 and 8 - 



Fig. 2. Effective Constellation at the relay for singular fade states, when the 
end nodes use BPSK constellation. 



To remove the distance shortening effect a procedure is 
given in [4-| when the nodes A and B use QPSK signal set. The 
procedure suggested in [4| to obtain the channel quantization 
and the clusterings, was using a computer algorithm, which 
involved varying the fade state values over the entire complex 
plane, i.e., < 7 < 00, < 9 < 2tt in small discrete steps 
and finding the clustering for each value of channel realization. 
But such an approach have many issues. In [4|, it is claimed 
that the clustering used by the relay is indicated to A and B 
by using overhead bits. However, the procedure suggested in 
(4) to obtain the set of all clusterings, was using a computer 
search algorithm (called Closest Neighbour Clustering (CNC) 
algorithm), which involved varying the fade state values over 
the entire complex plane, i.e., < 7 < 00, < 8 < 2ir and 
finding the clustering for each value of channel realization as 
discussed in previous sections. The total number of network 
codes which would result is known only after the algorithm 
is run for all possible realizations je^ 8 which is uncountably 
infinite and hence the number of overhead bits required is not 
known beforehand. Moreover, performing such an exhaustive 



search is extremely difficult in practice, especially when the 
cardinality of the signal set M is large. 

The implementation complexity of CNC suggested in O 
is extremely high: It appears that, for each realization of the 
singular fade state, the CNC algorithm of |]4) needs to be run 
at R to find the clustering. 

In the CNC algorithm suggested in [4], the network coding 
map is obtained by considering the entire distance profile. The 
disadvantages of such an approach are two-fold. 

• Considering the entire distance profile, instead of the 
minimum cluster distance alone which contributes dom- 
inantly to the error probability, results in an extremely 
large number of network coding maps. For example, for 
16-QAM signal set, the CNC algorithm results in more 
than 18,000 maps 0. 

• The CNC algorithm tries to optimize the entire distance 
profile, even after clustering signal points which con- 
tribute the minimum distance. As a result, for several 
channel conditions, the number of clusters in the clus- 
tering obtained is greater than the number of clusters in 
the clustering obtained by taking the minimum distance 
alone into consideration. This results in a degradation in 
performance during the BC phase, since the relay uses a 
signal set with cardinality equal to the number of clusters. 
For example, for 16-QAM signal set, the relay has to use 
signal sets of cardinality 16 to 29 flU. 

In |4), to overcome the two problems mentioned above, 
another algorithm is proposed, in which for a given je 1 , 
an exhaustive search is performed among all the network 
coding maps obtained using the closest-neighbour clustering 
algorithm and a map with minimum number of clusters is 
chosen. The difficulties associated with the implementation of 
the CNC algorithm carry over to the implementation of this 
algorithm as well. 

The contributions and organization of the paper are as 
follows: 

• A procedure to obtain the number of singular fade states 
for PAM and QAM signal sets is presented. 

• It is shown that for the same number of signal points M, 
the number of singular fade states for square A/-QAM 
is lesser in comparison with the number of singular fade 
states for M-PSK. The advantages of this result are two 
fold - QAM offers better distance performance in MA 
Phase and QAM requires lesser number of Latin squares 
(i.e., a reduction in number of overhead bits). 

• To remove the singular fade state (7 = 1,8 = 0) for 
VA/-PAM, a Latin Square is constructed. It is shown that 
the bit-wise XOR mapping cannot remove the singular 
fade state (7 = 1,0 = 0) for any M-QAM and a different 
mapping is obtained to remove the singular fade state 
(7 = 1,9 = 0), from the Latin Square to remove the 
singular fade state (7 = 1,0 = 0) for \/a7-PAM. 

• By simulation it is shown that the choice of 16-QAM 
leads to better performance for both the Rayleigh and 
the Rician fading scenario, compared to 16-PSK. 



The remaining content is organized as follows: 
In Section HI] we discuss the relationship between singular 
fade states and difference constellation of the signal sets used 
by the end nodes. We present expressions to get the number of 
singular fade states for PAM and square QAM signal sets in 
Subsections III- Al and lH-Bl respectivelv. In Subsection lH-Cl it is 
proved that the number of singular fade states for M -QAM is 
always lesser in comparison with that of A/-PSK signal sets. In 
Section [III] the clustering for a singular fade state is obtained 
through completing a Latin Square and a Latin Square for 
removing the singular fade state z = 1 is analytically obtained 
for PAM and QAM signal sets. In SectionHVl simulation results 
are provided to show the advantage of Latin Square scheme 
for QAM over XOR network coding scheme as well as Latin 
Square scheme for PSK signal sets under Rayleigh and Rician 
fading channel assumptions. 

II. Singular Fade states and Difference 
Constellations 

In this section we show the relationships between singular 
fade states and difference constellation of the signal set used 
by the end nodes. The following lemma discusses the location 
of singular fade states in complex plane for any constellation 
used at end nodes. 

Lemma 1: Let node A use a constellation S\ of size M\ and 
let node B use a constellation S2 of size M^. Let xa, x' a £ S\ 
and x B ,x' B £ S2, then the singular fade states z = je^ 6 are 
given by 



ij 



1 



je 



j<> 



XA 



(5) 



x B - x B 

where xa,x' a £ Si and x B ,x' B £ S2. 

Proof: The pair (xa,x' a ) and (x B ,x' B ) result in the same 
point in the effective constellation at the relay if the complex 
numbers xa + J&* x B and x' A + ^e^x^ are the same. The 
expression (0 is obtained by equating these complex numbers. 

■ 

From Lemma[T]it can be seen that all the singular fade states 
in the complex plane is of the form of ratio of difference 
constellation points of the signal sets used by end nodes, 
i.e., the singular fade states are decided by the difference 
constellation points. Henceforth, throughout the paper, we 
assume both the end nodes use same constellation, S. Let 
AS denote the difference constellation of the signal set used 
at the end nodes S, i.e., AS — {x{ — a^|xj,a^ £ S}. For a 
fade state z — je^ 6 to become a singular fade state, it has 
to satisfy (f5]), or in other words z(x' B — x B ) = (xa ~ x' A ), 
where (x' B —x B ) and (xa — x' A ) are any point in AS. Hence, 
a singular fade state can be alternatively defined as follows. 

Definition 3: A singular fade state z is a mapping Z from 
AS to the complex plane C so that at least one dk £ AS is 
mapped to some di £ AS. The set of all singular fade states 
is given by {Z : AS -> C|3 Z(d k ) = d t }. 

Remark 1: The singular fade state z = 1 is the mapping 
from AS to itself that maps every point to itself. 

In the rest of this section, we focus on the singular fade 
states for symmetric PAM and square QAM signal sets. 



S = -(vH-l)+2n -3-11 3 

(a) \/~M PAM constellation 



AS = -2(v / A7-l)+2n -6 



(b) Difference Constellation 



(c) Singular fade states 



Fig. 3. yM PAM constellation, difference constellation and singular fade 
states for \f~M = 4 



A. Singular Fade States of PAM signal sets 

The symmetric \/M-PAM signal set is given by 

S = -(VM-l)+2n, n£ (0, • • • ,VM-1) 

and its difference constellation is given by 

AS = -2(Vm - 1) + 2n, ue(0,- ,2(\/M-l)). 

For example the 4-PAM signal set and it's difference constel- 
lation is given in Fig |3(a)| and Fig |3(b)| respectively. For each 
of the difference constellation point, the pair in the signal set 
which correspond to this point is also shown. We will often 
consider only the first quadrant of A5 only, denoted as AS + , 
which for a general complex signal set is given by 

AS + = {a : real(a) > 0, imaginary(a) > 0}. 

The following lemma gives the number of singular fade states 
for PAM signal sets. 

Lemma 2: The number of singular fade states, for a regular 
\/A7-PAM signal set, denoted by N,^j_ pAM , is given by 



^(Vm-pam) — ^ ' 



Vm-i 

E 



p\n 



1 

1- - 

P 



(6) 



where p\n stands for prime number p dividing n. 

Proof: There are 2(yM — 1) non-zero signal points in 
the difference constellation AS and since AS is symmetric 
about zero there are \[M — 1 signal points in A5 f+ . All these 
are scaled version of nonzero elements of Z^g. 



TABLE I 

Singular fade states for 8-PAM 



n, Elements 
in AS+ 


Relative primes 
less than n 


Singular fade 
states, z > 1 


V>(n) 


1 









2 


1 


2 


1 


3 


U 


34 


2 


4 


1,3 


4 


2 


5 


1,2,3,4 


c b b 5 
J - 2 ' 5 ' 4 


4 


6 


1,5 


ft 


2 


7 


1,2,3,4,5,6 


7 7 7 7 7 7 
'-2> 3> 4' K> fi 


6 



The number of positive integers less than or equal to n that 
are relatively prime to n is given by Euler's totient function, 



if)(n) 



p\n 



1 



where the product is taken over distinct prime numbers p 
dividing n. To get the total number of relatively prime pairs in 
Z /jj, we take the sum over all nonzero n € ^JJS which gives 
Y.n^i' 1 n Y\ P \n {}- |) ■ ° ne relatively prime pair (a,b) 
gives two singular fade states, a/b and b/a. The multiplication 
factor 4 in © accounts for the negative side of the in-phase 
axis as well as the inverses. Finally, 2 is added to count the 
singular fade state z = 1 and z = — 1, ■ 

Example 2: Consider the case of 4-PAM (M = 16) signal 
set as given in Fig|3] There are 2(vM — 1) = 6 non-zero 
signal points in the difference constellation. Scaled A5 + is 
having (\/A7 — 1) = 3 signal points-{l, 2, 3}. And there are 
14 singular fade states- 

112 3 -1-1-2 -3 

'2'3'3' ' '2'~ '"2"'"3"'1~'~ ' '~2~ 



These singular fade states are shown in Fig j3(c)| Calculating 
Euler totient function ip(n) for n — 1,2,3 we get 0,1,2 
respectively and substituting in ©, leads to N^-pam) = 
2 + 4(0 + 1 + 2) = 14. 

Example 3: For 8-PAM signal set the singular fade states 
with z > 1 are shown in Table]]] For each such z given in the 
table there exists singular fade states — z , - and — -. Hence, 
totally, there are 70 (2 + 4(0 + 1 + 2 + 2 + 4 + 2 + 6)) singular 
fade states. 

B. Singular Fade States for QAM signal sets 

We consider square Af-QAM signal set S = {A m i + 

jA m Q} where A m j and A m Q take values from the v^f-PAM 

signal set -(y/M-l) + 2n, n € (0, • • • ,v / M-l). We use 
the mapping p, : S — > Z m given by 



A mI +jA mQ -> -[(VM-l+A mI )VM+(VM-l+A mQ )] (7) 

for concreteness and our analysis and results hold for any 
map. The difference constellation AS of square QAM signal 
sets form a part of scaled integer lattice with (2yM — l) 2 
points. The 16-QAM signal set with the above mapping and 
its difference constellation is shown in Fig |4(a)| and in Fig |4(b)| 



3 • 


7 • 


11 • 


15. 


2 • 


6 • 


10. 


W • 










1 • 


5 • 


9 • 


13. 


• 


i • 


8 • 


12. 



(a) 16— QAM constellation 



(b) The Difference Constellation 
Fig. 4. 16— QAM constellation and its difference constellation 

Definition 4: lfl2l The Gaussian integers are the elements 
of the set Z[j] = {a + bj : a, b € Z} where Z denotes the set 
of integers. 

The signal points in the difference constellation are Gaus- 
sian integers. To get the number of singular fade states for 
square QAM signal sets, the notion of primes and relatively 
primes in the set of Gaussian integers is useful. 

Definition 5: lfl2l A Gaussian integer a is called a Gaus- 
sian prime if and only if the only Gaussian integers that divide 
a are: 1, — l,j, —j, a, —a, aj and — aj. The Gaussian integers 
which are invertible in Z[j] are called units in Z[j] and they 
are ±1 and ±j. Let a, /3 € Z[j]. If the only common divisors 
of a and j3 are units, we say a and f3 are relatively prime. 

Lemma 3: The number of singular fade states for the square 
M-QAM signal set, denoted by Nm-qam is given by 



N M - 



M-QAM 



= 4 + 8(j)(AS + ) 



where 4>(AS + ) is the number of relative prime pairs in AS + . 



o. 



. 2 



(a) 4— QAM constellation 

i A5+. 



TABLE II 

Prime factors of Gaussian integers in AS+ 



Elements in AS+ 


Prime factors 


No. of relatively prime pairs 


1 


1 


11 


1+j 


1+j 


6 


2 


1+j 


6 


1+2} 


l+2j 


10 


2+j 


2+j 


10 


2+2j 


1+j 


6 


3 


3 


10 


3+j 


(l+j),(l+2j) 


5 


l+3j 


(l+j),(2+j) 


5 


3+2j 


3+2j 


11 


2+3j 


2+3j 


11 


3+3j 


(l+j),3 


5 





(b) The Difference Constellation 
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(c) Singular fade states 

Fig. 5. 4— QAM constellation, its difference constellation and singular fade 
states 



Proof: All the possible ratios of elements from AS give 
singular fade states. We consider only ratios in AS + and 
multiply the number of possible such ratios with a factor of 4 
to account the ratios with points in all the other quadrants. To 
avoid multiplicity while counting we take only relative prime 
pairs in AS + and one such pair (a, 6) gives two singular fade 
states a/b and b/a. Because of this the multiplication factor 
becomes 8. Finally, the factor 4 is added to count the units. 



Example 4: For a 4-QAM signal set shown in Fig |5(a)| the 
number of singular fade states, N(z4-qam) is given by 4 + 
8.1 = 12. Scaled AS + have only two elements {1, 1 + j} in 
this case as shown in Fig |5(b)| They form one relatively prime 
pair. The singular fade states ±1, ±j, ±1 ± j, , * . is shown 
in Fig. |5(c)| 

Example 5: Consider the case of 16-QAM signal set. Table 
[H] discusses the prime factorization of the elements in AS + . 
From the table there are 96 relatively prime pairs, but it counts 



the pair (a, b) and (b, a) separately. So there are 48 distinct 
pairs of relative primes, and from Lemma [3] Niq-qam turns 
to be 388. The singular fade states of 16-QAM is shown in 
FigH 

TABLE III 
Comparison between M-PSK and M-QAM on number of 

SINGULAR FADE STATES 



M 


No. of singular fade 
states for M-PSK 


No. of singular fade 
states for M-QAM 


4 


12 


12 


16 


912 


388 


64 


63,552 


8388 



C. Singular fade states of M-PSK and M-QAM signal sets 

In this section we show that the number of singular fade 
states for M-QAM signal sets is lesser in comparison with 
that of M-PSK signal sets. The advantages of this are two 
fold- QAM offers better distance performance and it requires 
lesser number of overhead bits since the required number of 
relay clusterings are lesser in the case of QAM compared with 
PSK. 

Lemma 4: The number of singular fade states for M- 
QAM signal set is upper bounded by 4(n 2 — n + 1), where 

n = -, which is same as 4(M 2 ) - (2M - 

i)Vm + i). 

Proof: There are [(2vM — l) 2 — 1] non zero signal points 
in AS which are distributed equally in each quadrant, i.e., the 

number of signal points in A S + , which 

we denote by n. The maximum number of relatively prime 



pairs in a set of n Gaussian integers is 



i(ra-l) 



Since an upper 



bound is of interest we substitute this in Lemma [3] instead of 
4>(AS + ). This completes the proof. ■ 

The number of singular fade states for M -QAM signal set is 
lesser in comparison with that of M-PSK signal sets. In fTTtl 
it is shown that the number of singular fade states for M- 
PSK signal set is M(^ - ^ + 1), in C(M 3 ). From Lemma 
|4] an upper bound on the number of singular fade states for 
M-QAM is in C(M 2 ). 

Example 6: The singular fade states of 16-PSK signal set 
is given in Fig |6(b)| There are 912 singular fade states in total. 
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(a) 16-QAM (b) 16-PSK 

Fig. 6. Singulai' Fade States for 16— QAM and 16— PSK modulation schemes 



The advantage of square QAM constellation is highly ef- 
fective in higher order constellations, for example 64-QAM 
is having 8,388 singular fade states where as a 64-PSK has 
63,552 singular fade states and relay has to adaptively use 
63,552 clusterings. With the use of square QAM constellations 
the complexity is enormously reduced. 

Ill, Exclusive Law and Latin Squares 

Definition 6: |7) A Latin Square L of order M with the 
symbols from the set Z t = {0, 1, • • • , t — 1} is an M x M 
array, in which each cell contains one symbol and each symbol 
occurs at most once in each row and column. 

In iflOl it is shown that when the end nodes use signal sets 
of same size all the relay clusterings which satisfy exclusive- 
law can be equivalently representable by Latin Squares, with 
the rows (columns) indexed by the constellation point used by 
node A (B) and the clusterings are obtained by taking all the 
slots in Latin Squares which are mapped to the same symbol 
in one cluster. 

A. Removing Singular fade states and Constrained Latin 
Squares 

The minimum size of the constellations needed in the BC 
phase is M, but it is observed that in some cases relay may 
not be able to remove the singular fade states with t = M and 
t > M results in severe performance degradation in the MA 
phase 0. Let (k, l)(k', I') be the pairs which give same point 
in the effective constellation Sr at the relay for a singular 
fade state, where k,k',l,l' € {0, 1, ...., M - 1} and k, k' 
are the constellation points used by node A and I, I' are the 
corresponding constellation points used by node B. If they are 
not clustered together, the minimum cluster distance will be 
zero. To avoid this, those pairs should be in same cluster. This 
requirement is termed as a singularity-removal constraint. So, 
we need to obtain Latin Squares which can remove singular 
fade states and with minimum value for t. Towards this end, 
initially we fill the slots in the M x M array such that for 
the slots corresponding to a singularity-removal constraint the 



same element is used to fill slots. This removes that particular 
singular fade state. Such a partially filled Latin Square is called 
a Constrained Partially Latin Square (CPLS). After this, to 
make this a Latin Square, we try to fill the other slots of the 
CPLS with minimum number of symbols. 

Definition 7: A Latin Square L T is said to be the Transpose 
of a Latin Square L, if L T (i,j) = L(j,i) for all i,j € 
{0,1,2,..,M-1}. 

Lemma 5: For any constellation, if the Latin Square L 
removes the singular fade state z then the Latin Square L T 
will remove the singular fade state z~ l . 

Proof: Let the singular fade state z as given in ([5]) with 
constraint {(xa,xb), (x' a ,x' b )}. Then, by taking the inverse 



x B - x B 



xa x a 
Now the constraints are modified to {(xb,xa),(x' b ,x' a )}, 
i.e., the role of node A and node B are interchanged, which 
clearly results in the transpose of the Latin Square. ■ 

From the above lemma, it is clear that we have to get Latin 
Squares only for singular fade states \z\ < 1 or \z\ > 1. 

The square QAM signal set has a symmetry which is ir/2 
degrees of rotation. This results in a reduction of the number of 
required Latin Squares by a factor 4 as shown in the following 
lemma. 

Lemma 6: If L is a Latin Square that removes a singular 
fade state z, then there exist a column permutation of L such 
that the permuted Latin Square L' removes the singular fade 
state ze j * /2 . 

Proof: For the singular fade state z is given in (0 with 
constraint {(xa, xb), (x' a , x b )}, the singular fade state ze^l 2 
is given by 



ze 



jir/2 



[XA 



ze 



j-k/2 



[x' B - x B ] 
[xa - x' A ] 



Al_ e Jir/2 



x n e 



-J't/2 



XB& 



-Jir/2] 



Since in the square QAM constellation there exist signal points 
with x'ge"^/ 2 and iBe^ 1 / 2 , all the constraints are changed 
but the new constraints are obtainable from the permutation 
of signal points in the constellation used by node B. The 
columns of the Latin Squares are indexed by the signal points 
used by B and the effective permutation in the constellation is 
representable by column permutation in the Latin Square. ■ 

Note that the fade state z — 1 or (7 = 1, = 0) is a singular 
fade state for any signal set. 

Definition 8: A Latin Square which removes the singular 
fade state z — 1 for a signal set is said to be a standard Latin 
Square for that signal set. 

When the signal sets is a 2 A -PSK signal set then, in 1 1 1 ] it 
has been shown that the Latin Square obtained by Exclusive- 
OR (XOR) is a standard Latin Square for any integer A. It turns 
out that for M— QAM signal sets the Latin Square given by 
bitwise Exclusive-Or (XOR) is not a standard Latin Square 
for any M > 4. This can be easily seen as follows: Any 
square M-QAM signal set (M > 4) has points of the form 
a, a + jb, a — jb, for some integers a and b. For z = 1, the 
effective constellation at R during the MA phase contains the 
point 2a can result in at least two different ways, since 2a = 
a + za = (a + jb) + z(a — jb) for z = 1. Let l\, l 2 and Is 
denote the labels for a, a + jb, and a ~ jb respectively. For 
the singular fade state z = 1, we have {(h,l\), (h>h)} as a 
singularity removal constraint. But the Latin Square obtained 
by bitwise XOR mapping does not satisfy this constraint since 



B. Standard Latin Square for \JM~-PAM 

In this subsection, we obtain standard Latin Squares for 
a/M-PAM signal sets. 

Definition 9: An M x M Latin square in which each row 
is obtained by a left cyclic shift of the previous row is called 
a left-cyclic Latin Square. 

Lemma 7: For a %/M-PAM signal set a left-cyclic Latin 
Square removes the singular fade state z = 1. 

Proof: Consider the \/M-PAM signal set with the signal 
points labelled from left to right as discussed in Section UU 
Let {(ki,lx)(k2,l2)} be a singularity removal constraint. To 
get the same point in the received constellation at the relay R, 
when z = 1, we have fci+Zi = k 2 + l 2 . Consider the following 
two cases satisfying this equality. Case (i): k 2 = h,l 2 = &i 
In this case the constraint becomes {(ki,li)(h,ki)}, i.e., the 
Latin Square which removes z = 1 should be symmetric about 
main diagonal. 

Case (ii): k 2 = fei + m, l 2 = h — m for any m < vM., 
The constraint now becomes {(fci, h)(k± + m,li — m)} which 
means the symbol in fci-th row and ^i-th column should be 
repeated in the k\ + 1-th row and the ^1 — 1-th column. 
It is easily seen that a left-cyclic Latin Square satisfies both 
this requirements. ■ 

Example 7: Consider the received constellation at the relay 
when the end nodes use 4-PAM constellation and let the 
channel condition be z = 1 as given in FigJT] The singularity 



removal constraints are 

{(0,1)(1,0)}, {(0,2)(1,1)(2,0)}, {(0,3)(1,2)(2,1)(3,0)}, 
{(1,3)(2,2)(3,1)}, and{(2,3)(3,2)}. 

The Latin Square which removes this singular fade state is 
given in Fig|9] 

(0,1X1,0) (0,3}(1,2)(2,1)(3,0) (2,3)(3,2) 
• • • • • » • 

(0,0) (0,2)(1,1)(2,0) (1,3)(2,2)(3,1) (3,3) 

Fig. 7. Received Constellation at the relay for 2 = 1. 
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Fig. 9. Left-cyclic Latin Square to remove the singular fade state z = 1 

C. Standard Latin Square for M—QAM 

In this subsection standard Latin Square for a square M- 
QAM constellation is obtained from that of \/M-PAM con- 
stellation. 

Let PAM - i, for i = 1, 2, • • • , \/M, denote the symbol 
set consisting of \[M symbols {(i - 1)VM, ((i - 1)VM) + 
1, ((i - \)7m) + 2, • • • , ((» - 1)VM) + (VM - 1)}. Let 
LpAM-i denote the standard Latin Square with symbol set 
PAM - i for VM-PAM and also let Lqam denote the 
standard Latin Square for A/-QAM. Then, Lqam is given 
in terms of Lpam-i, i = 1, 2, • • • , vAf , as the block left- 
cyclic Latin Square shown in Fig. [10] This is formally shown 
in the following Lemma. 

Lemma 8: Let PAM — i for i — 1, 2. • • ■ , \^M, denote 
the symbol set consisting of \f~M symbols {(? — 1)^/M , ((i — 
l)VM) + l,{(i-l)VM) + 2,--- ,((i-l)VM) + (>/M-l)} 
and let Lpam-i stand for the Latin Square that removes the 
singular fade state z — 1 with a symbol set PAM — i for 
%/M-PAM. Then arranging the cyclic Latin Squares Lpam-i 
as shown FigJTOl where each row is a blockwise left-cyclically 
shifted version of the previous row results in a Latin Square 
which removes the singular fade state 2=1 for M-QAM. 

Proof: Note that the matrix in Fig. [I0]is a Mx M matrix, 
which is also a \f~M~Y. \f~M block left-cyclic matrix where each 
block is a \[M x s/M left-cyclic matrix LpAM-i for some i. 

Let oi +jbi, a 2 +jb 2 ,a' 1 +jb[ and a 2 +jb' 2 , where a%, a\, bi 
and^e {-(\/M-l),-(\/M-3),--- ,(Vm~-3),(VM- 
1)} for i 6 {1,2} be four A/-QAM constellation points such 
that a\ + jb\ and a[ + jb[ are used by node A and a 2 + jb 2 
and a' 2 + jb' 2 are used by end node B, and result in a same 
point in the effective received constellation at the relay node 
for singular fade state z = 1, i.e., 



ai + jbi +a 2 + jb 2 = a[ + jb[ 



Jb' 2 - 



mi 



and b\ 



Let a[ = a\ 
{-2(VM-l),-2(VM-2)~ 



h 
,2( 



m 2 where mi , m 2 € 
M-2),2(VM-1)}. 
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Fig. 8. Standard Latin Square Lqam for 16-QAM. 



Then, a' 2 = a 2 — mi and 6 2 = ^2 ~ wi2- Then, using the map 
defined in [7] let 

fei = /i(ai + j'foi) 
/i = n{a 2 + jb 2 ) 

k 2 = n(a[ + jb[) = fi(ai + mi + j(6i + m 2 )) 
^2 = ^(4 + jfe 2 ) = /x(a2 - mj. + j{b 2 - m 2 )) 

Since, for z = 1, the four complex numbers result in 
the same point in the effective constellation at the relay, 
{(fei, h)(k 2 , l 2 )} is a singularity removal constraint for 2=1. 
From the above equations it follows that 

fa = fci + -(miVM + m 2 ) 
? 2 = Zi - -(miv / M + m 2 ) 

The above equations precisely mean the construction shown 
in Fig[l0] This completes the proof. 

■ 

The standard Latin Square for 16-QAM is shown in FigJH] 

We define a minimal Latin Square as, 

Definition 10: An M x M Latin Square with M symbols 
is termed as a minimal Latin Square. 

IV. SIMULATION RESULTS 

The proposed Latin Square (LS) Scheme ( IfTUl ) is based 
on removing the singular fade states. For 16-PSK all the 
912 singular fade states can be removed with minimal Latin 
Squares, but for 16-QAM some singular fade states cannot 
be removed with minimal Latin Squares. Since 16-QAM 
have only 388 singular fade states, in comparison with 912 
singular fade states of 16-PSK, 16-QAM offers better distance 
distribution in the MA stage. For a given average energy, the 
end to end BER is a function of distance distribution of the 
constellations used at the end nodes as well as at the relay. 



LpAM-2 



J PAM-y/M 



J PAAI-VM-1 



Fig. 10. Construction of Lqam for z = 1 



The simulation results for the end to end BER as a function of 
SNR is presented in this section for different fading scenarios. 

Consider the case when Ha,Hb,H' a and H' B are dis- 
tributed according to Rayleigh distribution, with the variances 
of all the fading links are assumed to be dB. The end to 
end BER as a function of SNR in dB when the end nodes 
use 16-QAM signal sets as well as 16-PSK signal sets with 
same average energy is given in FigJTT] The end to end BER 
for XOR network code for 16-QAM is also given. It can be 
observed that the LS Scheme for 16-QAM outperforms LS 
Scheme for 16-PSK as well as XOR network code. 

Consider the case when Ha,Hb,H' a and H' B are dis- 
tributed according to Rician distribution, with the Rician factor 
of 5 dB and the variances of all the fading links are assumed 
to be dB. In Fig[l2]the end to end BER as a function of SNR 
in dB for LS scheme for 16-PSK, 16-QAM and XOR network 
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Fig. 11. SNR vs BER for different schemes when the end nodes use 16-QAM 
and 16-PSK for a Rayleigh fading scenario. 



Fig. 12. SNR vs BER for different schemes when the end nodes use 16-QAM 
and 16-PSK for a Rician fading scenario with Rician factor 5 dB. 



coding for 16-QAM is given. It is observed that the LS scheme 
gives large gain over the XOR network coding scheme. The 
LS scheme for QAM is better in end to end BER performance 
in comparison with the LS scheme for PSK. 

V. DISCUSSION 

In this paper, for the design of modulation schemes for 
the physical layer network-coded two way relaying scenario 
with the protocol which employs two phases: Multiple access 
(MA) Phase and Broadcast (BC) phase, with both end nodes 
use square QAM constellation is studied. We showed that 
there are many advantages of using square QAM constellation. 
With the help of the relation between exclusive law satisfying 
clusterings and Latin Squares we propose a method to remove 
the singular fade states. This relation is used to get all the 
maps to be used at the relay efficiently. We proposed a 
construction scheme to get the Latin Square for square QAM 
constellation from PAM constellation. Here we concentrated 
only on singular fade states and the clusterings to remove that 
with only the minimum cluster distance under consideration. 
We are not considering the entire distance profile as done in 
|4). Unlike in the case of 0, we could remove most of the 
singular fade states with standard Latin Square and its isotopes. 
We presented the simulation results showing the end to end 
BER performance when the end nodes use PSK constellation 
as well as QAM constellations. 
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